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ABSTRACT 

Wc have investigated the eclipse of the 23 ms pulsar, PSR J0737— 3039A, by its 2.8 s companion 
PSR J0737— 3039B in the recently discovered double pulsar system using data taken with the Green 
Bank Telescope at 820 MHz. We find that the pulsed flux density at eclipse is strongly modulated with 
the periodicity of the 2.8 s pulsar. The eclipse occurs earlier and is deeper at those rotational phases 
of B when its magnetic axis is aligned with the line of sight than at phases when its magnetic axis is 
at right angles to the line of sight. This is consistent with the eclipse of A being due to synchrotron 
absorption by the shock- heated plasma surrounding B, the asymmetry arising from the higher plasma 
densities expected in the B magnetosphere's polar cusps. 

Subject headings: pulsars: general — pulsars: individual (J0737— 3039A, J0737— 3039B) — radiation 
mechanisms: non-thermal — binaries: general 



B. In their early analysis of the exploratory observations 
performed at the 100-m Green Bank Telescope (GBT) 
in 2003 December and 2004 January, Kaspi et al. (2004) 
found no evidence for modulation of the eclipse properties 
with B's rotational phase. We have analyzed this publi- 
cally available dataset and, in this Letter, show that the 
eclipse properties of A are in fact strongly dependent upon 
the rotational phase of B. 

2. OBSERVATIONS AND ANALYSIS 

The double pulsar system J0737-3039 was observed at 
the GBT in 2003 December and 2004 January using re- 
ceivers centered at at 427, 820 and 1400 MHz. The 5-hr 
820-MHz observation discussed here was acquired with the 
GBT Spectrometer SPIGOT card with a sampling time of 
40.96 fis and 1024 synthesized frequency channels covering 
a 50-MHz bandwidth. For further details of the observa- 
tions and data acquisition system, see Ransom et al. (2004) 
and references therein. 

The data were dedispersed and folded using freely avail- 
able software tools (Lorimer 2001) assuming the nominal 
dispersion measure of 48.9 cm^'^ pc (Burgay et al. 2003) 
and using ephemerides for A and B from Lyne et al. (2004). 
We folded the data with 512 phase bins and, for each A 
pulse, calculated the mean pulsed flux density in two on- 
pulse phase windows centered on the two pulse compo- 
nents (Burgay et al. 2003). From these, we subtracted the 
baseline level calculated from two off-pulse windows. The 
on- and off-pulse regions contained 95 and 320 phase bins, 
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1. INTRODUCTION 

The two pulsars in the recently discovered 
double pulsar binary system, J0737— 3039A and 
J0737-3039B (hereafter simply "A" and "B") have pe- 
riods of Pa = 23 ms and Pb = 2.8 s. They are in a 
2.4-hr mildly-eccentric orbit that we view nearly edge-on, 
with an inclination angle of 87° ± 3° (Burgay et al. 2003, 
Lyne et al. 2004). The phenomenology exhibited by this 
system is extremely rich. The flux density of the B pulsar 
varies dramatically and systematically around the orbit 
(Lyne et al. 2004, Ramachandran et al. 2004, McLaugh- 
lin et al. 2004), indicating significant interaction between 
the two pulsars. While the flux density of the A pulsar 
is constant for most of the orbit, for ~ 30 s around its 
superior conjunction, A is eclipsed by the magnetosphere 
of B (Lyne et al. 2004). This eclipse is asymmetric, with 
the flux density of A decreasing more slowly on ingress 
than it increases on egress. The duration of the eclipse is 
mildly dependent on radio frequency, with eclipses lasting 
longer at lower frequencies (Kaspi et al. 2004). 

Theoretical models by Arons et al. (2004) and Lyutikov 
(2004) explain these eclipse properties in the context of 
synchrotron absorption of A's radio emission by the shock- 
heated plasma surrounding and containing B's magneto- 
sphere. Because B's magnetosphere will respond differ- 
ently to the pressure of the relativistic wind of A at differ- 
ent rotational phases of B, one prediction of these models is 
that the density of absorbing plasma, and hence the prop- 
erties of A's eclipse, will depend on the rotational phase of 
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respectively. For each orbit, we created a light curve of 
A"s pulsed flux density with time by averaging every f2 
pulses, for an effective time resolution of ~ 0.27 s. Orbital 
phases were calculated using the modified Julian dates of 
pulse arrival at the Solar-System barycenter and the pul- 
sar ephemeris. Analysis of these light curves reveals that, 
for the majority of the orbit, the pulsed flux density of A 
is constant within measurement uncertainties, exhibiting 
no obvious orbital phase-dependent variations. However, 
close to superior conjunction (i.e. orbital phase 90°, where 
orbital phase is defined as the longitude from the ascending 
node), the pulsed flux density of A varies dramatically. 
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Fig. 1. — The pulsed flux density of A versus time (with respect 
to superior conjunction) and orbital phase for (top three panels) the 
three eclipses in the 820-MHz observation and (bottom panel) all 
three eclipses summed. In the individual eclipse light curves, ev- 
ery 12 pulses have been averaged for an effective time resolution of 
~ 0.27 s. Every 100 pulses have been averaged to create the lower, 
composite light curve for an effective time resolution of ~ 2.3 s. 
Pulsed flux densities have been normalized such that the pre-eclipse 
average flux density is unity. 



Figure 1 shows the light curves of A for all three eclipses 
included in the 820-MHz observation. In the lower panel, 
we also show a composite light curve, averaged over all 
three eclipses and over 100 A pulses, for an effective time 
resolution of ^ 2.3 s. Our measurements of the eclipse 
duration and of ingress and egress shapes from this com- 
posite light curve are consistent with the results of Kaspi 
et al. (2004), who averaged data in 2 s intervals. They cal- 
culated an eclipse duration (defined by the full width at 
half maximum of the light curve) of 27 s. which was used 
to place a limit of 18,600 km on the size of the eclipsing 
region. They also found that eclipse ingress takes roughly 
four times longer than egress and that the eclipse duration 
is frequency dependent, lasting longer at lower frequencies. 
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Orbital phase (degrees) 

Fig. 2. — As in Figure 1, the top three panels show the light curves 
for the three individual eclipses, but covering a narrower range of or- 
bital phase. The vertical dashed lines indicate the measured arrival 
times of the pulses of B. The lower panel shows the light curve aver- 
aged over all three eclipses of A, with the second and third eclipses 
shifted by < itPs/2 so that the B pulses arrive at the same phase 
as during the first eclipse. Because limited signal-to-noise ratio pro- 
hibits us from summing over a smaller number of pulses, it is likely 
that the true structure of the flux density peaks and dips is finer 
than it appears in this figure. 

Inspection of the top three panels of Figure 1 shows that 
the rise and fall of the flux density during eclipse is not 
monotonic. In fact, up until orbital phases approaching 
90°, the light curves of A show amplitude peaks consis- 
tent with pre-eclipse levels. Figure 2 shows an expanded 
view of Figure 1 covering a smaller range of orbital phases 
centered on the eclipse. On this plot, we also indicate the 
measured barycentric arrival times of the pulses of the 2.8 s 
pulsar B. This demonstrates clearly that the pulsed flux 
density of A is modulated in synchronism with 0.5 Pb- In 
all three eclipses, we see negative dips in pulsed flux den- 
sity occurring first at a time about 0.5 Pb out of phase 
with the B pulses, and later also at a time in phase with 
the B pulses. In order to determine more sensitively how 
the flux density of A depends upon the rotational phase 
of B, we shifted the light curves of the second and third 
eclipses by up to ± 0.5 Pb so that the phases of the B 
pulses were identical to those of the first eclipse. We then 
summed the light curves of all three eclipses to create the 
bottom panel of Figure 2. Dividing each 2.8 s window of 
B's rotational phase into four equal regions, we can calcu- 
late an average light curves for each region, as shown in 
Figure 3. The light curves for each of the B pulse phase 
windows vary smoothly. 

Eclipse durations range from 20 s to 34 s for the four 
B phase windows, calculated, as in Kaspi et al. (2004), to 
be the full width at half maximum of the eclipse. Ingress 
occurs first at B phase 0.5 and then at 0.0, when, assum- 
ing B is an orthogonal rotator, the magnetic axis of B is 
aligned with the line of sight to A. It occurs later and 
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more gradually for the other two phase windows centered 
on phases 0.25 and 0.75, when the magnetic axis of B is 
at right angles to the line of sight. The orbital phases of 
egress are similar for all B pulse phases. For light curves 
at B phases centered on 0.0 and 0.5, the eclipses are es- 
sentially symmetric, with minimum beginning at orbital 
phases near 89.4° and ending at orbital phases near 90.6°. 
The light curves for B phases centered on 0.25 and 0.75 
are less symmetric, with the flux density minimum occur- 
ring briefly, around orbital phases 90.2° — 90.5°. Only at 
around orbital phase 90.5° is flux density of all phases of 
B consistent with being zero. 

It is important to recognize that the features we see in 
the 820-MHz data cannot be due to confusion with the 
single pulses of B (e.g. McLaughlin et al. 2004) since the 
pulsed features seen in Fig. 2 are not in phase with the 
radio pulses of B, and, in addition, the latter, which have 
width of ^ 50 ms, would be removed by the on-pulse/off- 
pulse subtraction procedure used to determine A's flux 
density. 




Orbital phase (degrees) 

Fig. 3. — Averaged light curves for four regions of B pulse phase, 
with a smooth curve calculated with a spline fitting procedure drawn 
through the individual flux density measurements. The four regions 
are centered on pulse phases 0.0, 0.25, 0.50 and 0.75, with, for ex- 
ample, the curve for 0.0 covering the 0.25 of pulse phase centered on 
B's radio pulse. Figure 4 shows the geometry of the system at this 
B pulse phase. The vertical bar indicates the typical ± one-sigma 
measurement error in each flux density. 



(6 X lO'^'^ erg s^^) is ^ 3600 times greater than that of 
the slower B pulsar. Because the dynamic pressure of A's 
wind is balanced against B's magnetic pressure only at a 
point well inside the light cylinder of B, A's wind con- 
fines the B magnetosphere on the side facing A, with a 
long, magnetospheric tail extending behind B. The colli- 
sion of A's wind with the magnetosphere of B leads to the 
formation of a bow shock consisting of hot, magnetized 
plasma. This plasma will surround the cometary shaped 
magnetosphere of B, creating what Arons et al. (2004) call 
a "magnetosheath" . In Figure 4, we show a schematic of 
the interaction of A's relativistic wind with B's magne- 
tosphere. The rotation of B inside the magnetosheath is 
expected to modulate the shape of the sheath and to pro- 
duce cusps within it, leading to higher plasma densities 
in the cusps. Note that the plasma will not reach all the 
way to the polar caps of B, but is confined at the "magne- 
topause" , where the plasma pressure equals the magnetic 
pressure. An excellent description of the shape of the mag- 
netosphere and of the polar cusps, is given in Arons & Lea 
(1976). 
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Fig. 4. — Cartoon (not to scale) showing the interaction between 
the relativistic wind of A and the magnetosphere of B when the radio 
beam of B is pointing towards the Earth (i.e. pulse phase 0.0). For 
a more detailed picture of this interaction, see Figure 2 of Arons et 
al. (2004). 



3. DISCUSSION 

Our analysis of these data differs from that of Kaspi 
et al. (2004) in one significant respect. While Kaspi et 
al. (2004) calculated light curves with 2 s time resolution, 
we have integrated over only 12 A pulses, affording us ^ 
7 times better time resolution. Kaspi et al. (2004) did 
detect some emission from A on short time scales during 
eclipse, but their 2 s integration time was not sufficient to 
detect the B modulation. We have shown that the eclipse 
duration varies considerably with B phase and, hence, that 
the properties of the occulting medium are in fact different 
in different regions of B's magnetosphere. 

The observed eclipse phenomenology can be understood 
in the context of a model invoking synchrotron absorp- 
tion of the radio emission from A in a magnetosheath sur- 
rounding B's rotating magnetosphere (Arons et al. 2004, 
Lyutikov 2004). Much of the rotational energy of the 
faster-spinning A pulsar is emitted in the form of a rela- 
tivistic, magnetized wind composed primarily of electron- 
positron pairs. The rotational energy loss rate of A 



As already discussed qualitatively by Kaspi et al. (2004) 
and quantitatively by Arons et al. (2004) and Lyutikov 
(2004), synchrotron absorption in such a magnetosheath 
can explain many of the properties of the A eclipse, includ- 
ing the sharp edge and the mild frequency dependence. If 
B is rotating prograde to the orbital velocity, it can also 
explain the asymmetry between eclipse ingress and egress. 
The observed changes in eclipse properties with B pulse 
phase give further support to the model. The pulses from 
A are first absorbed at the two B pulse phases (i.e. 0.0 
and 0.5) at which the magnetic axis of B is aligned with 
the line of sight to A. This can be attributed to the higher 
density of absorbing plasma in B's polar cusps. Note that 
the pulses at B pulse phase 0.50 (i.e. when the radio-bright 
pole is pointed towards A) are absorbed before those at B 
pulse phase 0.0. When the line of sight to A is at a right 
angle to B's magnetic axis, the pulses of A travel through 
regions of higher plasma density later in the orbit and start 
fading closer to superior conjunction, where the opacity of 
the magnetosheath tail could also become predominant. 
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As shown in Figure 3, the ingress phase and echpsc du- 
ration for these two B phases (i.e. 0.25 and 0.75) differ, 
with absorption happening later at phase 0.75. This may 
be due to the rotationally-induced asymmetry predicted 
by Arons et al. (2004). 

The above discussion, and the simulations of Arons et 
al. (2004), assume that B is an orthogonal rotator. As 
Demorest et al. (2004) and Arons et al. (2004) discuss, it 
is quite likely that this is the case as the wind torque from 
A should have aligned B's rotation axis with the orbital 
angular momentum. Furthermore, it is very difhcult to 
explain the similarity between the light curves at B pulse 
phases 0.0 and 0.5 and phases 0.25 and 0.75 if B is not 
nearly orthogonal. 

From the measured eclipse durations, it is clear that 
the size of the eclipsing region is much smaller than the 
light cylinder of B. However, because the plasma den- 
sity varies greatly depending on the rotational phase of 
B, determining the extent and density of the eclipsing re- 
gion is not trivial and will require detailed modeling of 
the plasma density variations within B's magnetosphere. 
Due to the effects of relativistic periastron advance and 
geodetic precession, the durations and morphologies of 
the eclipse and their dependence on B rotational phase 
shoiild change significantly with time. Forthcoming obser- 
vations will allow us to sensitively probe the dynamics of 
the wind/magnetosphere interaction and of the geometry 
of the system. These observations, at multiple frequen- 
cies, will also allow us to test the hypothesis of Arons et 
al. (2004) that there should be no eclipse at frequencies 
above 5 GHz, already supported by the measurements of 
Kaspi et al. (2004). 

Because pulsars emit a large portion of their spin-down 
energy in the form of a relativistic wind, understanding 
these winds is crucial for forming a complete picture of the 



pulsar energy budget. However, because the winds cannot 
be directly observed, progress in this field has been limited 
to studies of pulsar wind nebulae and bow shocks. If the 
eclipses of A are indeed caused by synchrotron absorption 
in B's magnetosheath, both Arons et al. (2004) and Lyu- 
tikov (2004) reach the very interesting conclusion that the 
density of A's wind must be at least four orders of mag- 
nitude greater than is expected given currently accepted 
models of pair creation. The Arons et al. (2004) model 
also implies a very low wind magnetization. However, in 
our paper describing the drifting features observed in the 
single pulses of B (McLaughlin et al. 2004), we concluded 
that that most of the spin-down energy is carried by the 
Poynting flux rather than by energetic particles. Clearly, 
self-consistent models which can explain both of these phe- 
nomena are needed. The rich phenomenology offered by 
the double pulsar system will offer us new ways to explore 
the physics of these winds on much smaller scales than has 
been possible with other systems. Eclipse measurements 
of the double pulsar system may also be important for 
constraining pulsar pair creation scenarios. 
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